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Plane symmetric self-similar solutions to Einstein's four-dimensional theory of gravity are studied 
and all such solutions are given analytically in closed form. The local and global properties of 
these solutions are investigated and it is shown that some of the solutions can be interpreted as 
representing gravitational collapse of the scalar field. During the collapse, trapped surfaces are 
never developed. As a result, no black hole is formed. Although the collapse always ends with 
■ " i" ' spacetime singularities, it is found that these singularities are spacelike and not naked. 
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^ . I. INTRODUCTION 

^~i , 

^^ ' Recently, self-similar solutions to the Einstein field equations have attracted lots of attention, not only because 

I [ the corresponding problem can be considerably simplified and, as a result, can be studied analytically [1], but also 

J> ■ because their relevance to critical phenomena in gravitational collapse, which were first discovered by Choptuik in 

l/~} ' 1993 in his numerical study of collapsing spherically symmetric scalar field [2]. The phenomena have lately been found 

04 in various matter fields collapse [3,4]. In particular, one of the present authors studied analytically the gravitational 

^D • collapse of cylindrically symmetric scalar field in four-dimensional spacetimes, and found a class of exact solutions with 

self-similarity [5]. It was shown explicitly that one of the solutions has precisely one unstable mode. By definition, 

this is a critical solution that sits on the boundary separating two different phases in the phase space of the initial 

(-^ [ data. This serves as the second analytical model of critical collapse found so far. The first analytical model was 

■»-> ' found by Garfinkle in 2001 in 2 -I- 1 gravity [6], right after the numerical simulations of Pretorius and Choptuik [7] 

^:Li and Husain and Olivier [8]. The linear perturbations of Garfinkle's self-similar solutions, denoted by S[n\, where n 

I is a parameter, were systematically studied in [9] and [10]. Due to different boundary conditions imposed, different 

wp." results were obtained. In particular, in [9] it was found that the solution with n = 3 has precisely one unstable mode, 

. . , while in [10] it was found that the one with n = 4 has one and only one unstable mode. Although the solution with 

^ • n = 4 is best matched with the numerical critical solution found by Pretorius and Choptuik, the exponent, 7, of the 

black hole mass 



o 



>< 



c5 ■ Mbh^{p-P*)\ (1-1) 



obtained by [11] 



'"W <'-^' 



is quite different from the numerical one, 7 ^ 1.2 [10], where ki is the unstable mode of the critical solution. On 
the other hand, the critical solution found in [9] is different from the numerical one [7], but the resultant exponent, 
7 — 4/3, is close to the numerical one found by Pretorius and Choptuik, although it is still quite different from the 
numerical one found by Husain and Olivier, 7 ~ 0.81 [8]. 

In this paper, we have no tendency to resolve the above disputations, but look for more analytical solutions that 
might represent critical collapse. In the analytical studies of critical collapse, the investigation is usually divided into 
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two steps: One first finds some particular solutions by imposing certain symmetries, such as, homothetic self-similarity 
(HSS). This can mathematically simplify the problem considerably. For example, in the spherically symmetric case, 
by imposing HSS symmetry the Einstein field equations will be reduced from PDE's to ODE's. Once the particular 
solutions are known, one can study their linear perturbations and find the spectra of the corresponding eigen-modes. 
If a solution has precisely one unstable mode, it may represent a critical solution, sitting on a boundary that separates 
two different basins of attraction in the phase space. In this paper, we shall restrict ourselves only to the first step. 
In particular, we shall study self-similar solutions of massless scalar field with plane symmetry in Einstein's four- 
dimensional theory of gravity, and present all such solutions in closed form. Then, we shall study their local and 
global properties, whereby show some of them can be interpreted as representing gravitational collapse of the scalar 
field. During the collapse, trapped surfaces are never developed. As a result, no black hole is formed. Although the 
collapse always ends with spacetimc singularities, it is found that these singularities are spacelike and not naked. 

II. PLANE SYMMETRIC SPACETIMES WITH SELF-SIMILARITY 

The general metric for spacetimes with plane symmetry can be cast in the form [12], 

ds^ ^ 2e-^^("^^)dudw - e-^("'") {dx^ + dy^) , (2.1) 

where {x'^} = {M,i',a;,y}, {^ = 0,1,2,3), and — cxo < x^ < +oo. The three Killing vectors that characterize the 
symmetry are given by ^(i) — dx, C(2) = dy and ^(3) — xdy — ydx- The non-vanishing components of the Ricci tensor 
are given by [13] 

(2.2) 

(2.3) 

(2.4) 

(2.5) 

where ( ).„ = d{ )/du and so on. A massless scalar field satisfies the Klein-Gordon equation, g°''^4';afi — 0, which in 
the present case takes the form, 

20,™ - (f>,uU^v - (t>,vU^u = 0. (2.6) 

However, this equation is not independent of the Einstein field equations and can be obtained from the Bianchi 
identities G^a;i3g°'^ = 0. 

Spacetimes with homothetic self- similarity (or self- similarity of the first kind) is usually defined by the existence of 
a conform Killing vector ^^ that satisfies the equations [14], 

where a semicolon ";" denotes the covariant derivative. It can be shown [15] that for the spacetimes described by 
Eq.(2.1) there are two types of self-similar solutions. The first type is given by 

M{u,v) ^ M{z), U{u,v)^ S{z)-2ln{-u), 

^f^d^^udu + vd,, z = -^, (2.8) 

i-u) 

while the second type is given by 

M{u, v) = M(z), U{u, v) = U{z), 

ed^^x^dx, z=j^. (2.9) 

Solutions with the properties of Eq.(2.8) will be referred to as the Type A self-similar solutions, and the ones with 
the properties of Eq.(2.9) will be referred to as the Type B self-similar solutions. In the following we shall consider 
them separately. For the sake of convenience, we shall first restrict ourselves to the region u < 0, w > 0. Once we 
find the self-similar solutions in this region, we shall extend them to other regions whenever it is necessary. 



^uu 


= - {2U,uu - U^J + 2U.uM,u) 
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Rvv 
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Rxx 


= Ryy = -e^^-^ {U,UV - U^uU, 



III. SOLUTIONS OF MASSLESS SCALAR FIELD WITH TYPE A SELF-SIMILARITY 

The solutions with Type A self-similarity are characterized by Eq.(2.8), for which the non-vanishing components of 
the Ricci tensor given by Eqs.(2.2)-(2.5) become 

(3.1) 
(3.2) 
(3.3) 

(3.4) 

where a prime denotes the ordinary differentiation with respect to z. For a massless scalar field, the Einstein field 
equations read 
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(3.5) 



where k[= SttG/c^] is the Einstein coupling constant. In this paper we shall choose units such that k — 1. It can be 
shown that a massless scalar field </)(«, v) that is consistent with spacetimcs with homothetic self-similarity must take 
the form. 



(u, v) = 2q In(-u) + f{z), 



(3.6) 



where q is an arbitrary constant, and (p{z) is an arbitrary function, which will be determined by the Einstein field 
equations (3.5). Inserting Eqs.(3.1)-(3.4) and (3.6) into Eq.(3.5), we find that 



z{S" -^S'^\ +(zS' + 2)M' 
I 



{Zip' - 2qf 



z [M" + S") - -zS' + M' = if' [zLp' - 2q) , 
S" - ^5'' + S'M' = ip'\ 

^-S'-'- = 0. 

b' z 



Eq.(3.10) has the general solution. 
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where a and b are two integration constants. When ah — 0, it can be shown that the corresponding solutions 
correspond to null dust fluid, and in this paper we shall not consider this case. Thus, from now on we shall assume 
that ah ^ 0. Substituting the solution (3.11) into Eqs.(3.7)-(3.9), we find that 



M' = 



1 



2az (az'^ + b) 
a {zifi' — 2q) + bip' = 

, ^// , ,1 2abz 
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{az^ + b) 
On the other hand, one can show that Eq.(2.6) takes the form 
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(3.12) 
(3.13) 
(3.14) 

(3.15) 



It can be shown that Eq.(3.14) is not independent. As a matter of fact, it can be obtained from Eqs.(3.12) and (3.15). 
Thus, in the following we need only to consider Eqs.(3.12), (3.13) and (3.15). Integrating Eq.(3.15) we find that 



^'(z) ^ 



2aqz + Ci 



az-^ 



where ci is a constant. Inserting Eq.(3.16) into Eq.(3.12), we find that the constant ci must be given by 
To have Ci real, we must require ab < 0. Then, from Eqs.(3.12) and (3.16) we obtain. 
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(3.18) 



where Mq and ipo are other integration constants, and 



±[2(l + 2g2)] 

5X1/2 

a 



.)]i/2^ 
, ab <0. 



(3.19) 



Rescaling the coordinates x and y, without lose of generality, we can always set |a| = 1. Then, the solutions can be 
finally written in the form. 



M{z) = 2q^\n 



S{z) = — In a — z 



— eg In 
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(f{z) = q\n\a — z \ -\ — cln 
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(3.20) 



To study the above solutions, let us first consider the region where a > z. Then, it can be shown that 

[a[—u) — v\ [a(—u) + v\ 

X [cq (v'^ + a^u^) - a{l + Aq^) uv\ , 



e*^ = 



(3.21) 



Since 2q^ ± eg + 2 is always positive, from the above expressions we can see that the spacctime is singular on the 
hypersurface v = a(—u). The normal vector to this surface is given by 



_ d(v + au) 
" = d^^ = '5m + "^m' 



(3.22) 



which is always timelike, as now we have n^n" = 2ae^ > 0. Therefore, the singularity is spacelike and services as 
the up boundary of the region v > 0, m < 0. On the other hand, from these expressions we can also see that the 
spacetime is free of curvature singularity on the hypersurface w = 0, although the metric coefficient M is singular. 
Thus, to have a geodesically maximal spacetime, the solutions need to be extended beyond this surface. For the sake 
of convenience, in the following we shall consider the two cases 2g2 < l and 2g2 > 1 separately. 



A. a>z, 2g^ < 1 
In this case, introducing two new coordinates u and v via the relations 

/ -\n -n 

U = — ( — u) , V — V , 



(3.23) 



where 



1 



■'- I- 2g2 
we find that in terms of u and v the metric takes the form, 



(3.24) 



ds^ 



[ai-uY 



a(-u)" + i;"]''^ -"Uudv 



{a^u^^ - z;2") {dx^ + dy^) , 



(3.25) 



from which we can see that the metric coefficients are no longer singular on the hypersurface v = or i; = 0. The 
coordinate transformations (3.23) map the region w > 0, u < into the region w > 0, u < 0, which will be referred 
to as Region /, as shown in Fig. 1. Then, the region t; < 0, u < represents an extended region, which will be 
referred to as Region //. It should be noted that this extension is not always physically acceptable. In particular, the 
extension may not be analytical and the extended metric coefficients may not be real, unless n is an integer. When 
the extension is not analytical, it is also not unique. Therefore, to have an unique extension, in the following we shall 
consider the case where n is an integer. Then, in terms of v and u we find that 



4ae**" {-uv) 



X [cq (f;2" + a2fl2») + a {l + 4q^) {-uvf] , 



(3.26) 



from which we find that (/)^q, is timelikc only when n is an odd integer in the extended region, //. The spacetime is 
always singular on the hypersurface v = a~^'"u. This singularity is timclike, and the corresponding Penrose diagram 
is given by Fig. 1. Due to this singular behavior, it is found difficult to interpreted the corresponding solutions as 
representing physically reasonable model of gravitational collapse. 




FIG. 1. The Penrose diagram for the solutions given by Eq. (3.25) with n(> 1) being an integer. Region / is defined as 
V > 0, M < 0, and Region 77 as w < 0, u < 0. The spacetime is singular on the both horizontal and vertical double lines ab 
and ac, where v = a{—u) and v — a^^'"'u, respectively. The line be represents the past null infinity u = — oo. 



B. a>z, 2g^ > 1 



In this case, the spacetime in the region w > 0, u < is already geodesically maximal, and docs not need to be 
extended beyond the hypersurface v ^ 0. To see this, let us consider the null geodesies along u, x, y — Const., which 
are given by 



d'^v dM{uo,v) f dv 



dA2 



dv 



dX^ =0' 



(3.27) 



where A denotes the afSne parameter along the null geodesies. Near the hypersurface u = 0, it can be shown that the 
above equation has the following general solution, 
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0C1A+C2 
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(3.28) 



where ci.2 are the integration constants. Thus, as w ^ 0, we have |A| — > cxo, that is, the "distance" between the point 
{u,v) — (mo,0) and any other one along the null geodesic u = uq, say, (u,w) = {uq,vq), where vq > 0, is infinitively 
large. In other words, the hypersurface w = in the present case represents a physical boundary of the spacetime. 
Then, the corresponding Penrose diagram is given by Fig. 2. 




FIG. 2. The Penrose diagram for the solutions given by Eq. (3.20) with a > z, 2q^ > 1. The spacetime is singular on the 
horizontal double line ab, but not on the ones ac and be that now serve as the null infinity boundaries of the spacetime, where 
the line ac is given by i" = and the one bchy u — —00. 

In the region m < 0, v > 0, Eq.(3.21) shows that the scalar field is timelike. Thus, one may consider the above 
solution as representing gravitational collapse of the scalar field, and the collapse always forms a spacetime singularity 
on the hypersurface v — a{—u). Although no trapped surface is formed, the singularity is not naked, as it can be 
seen from Fig. 2. 

IV. SOLUTIONS OF MASSLESS SCALAR FIELD WITH TYPE B SELF-SIMILARITY 



For the Type B self-similar solutions, we find that the non-vanishing components of the Ricci tensor are given by 

Run = ^\z (u" - lu'^) + (2 + zM') U'] , (4.1) 

i?„„ = ^ z[M" + U" -^U''^]+M' + U' , (4.2) 
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R.. - Ryy = -^e''-^ fu" - W' + -] . 



U^ 



z 



(4.3) 
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One can show that the corresponding scalar field has to take the same form, Eq.(3.6). Since (f) = 4'{u, v), we find that 
Rxx ~ <i>\ ~ 0; which together with Eq.(4.4) yields 



9 U' 



(4.5) 



This equation has the general solution 

t/=-ln[aln(z) + /3], (4.6) 

where a and (3 are two integration constants. On the other hand, it can be shown that now Eq.(2.6) becomes, 



6 



(z(p')' + (9 - zf') U' - 0. (4.7) 
The general solution of it is given by 

ip{z) = <7ln(z)+7ln[Q;ln(z)+/3] + ipo, (4.8) 

where 7 and ipo are arbitrary constants. Inserting Eqs.(3.6), (4.3), (4.6) and (4.6) into the Einstein field equations for 
the OT-component, we find that 

M' ^ '-^ + "/;:^7\ - ^ [.ln(z) + P] , (4.9) 

z 2z[a\ii{z) + 13] az ^ ^ -^ ' ^ ^ 



which has the solution 



^2 



M(z) = (I - 2q7)ln(z) + i (f - 272) ln{[aln(z) +/3]} - ^ [a\n{z) + pf + Mq, (4.10) 

where Mq is another integration constant. Substituting the solutions into the rest of the Einstein field equations, we 
find that the constant 7 must be given by, 

7=1 (4.11) 

Then, the general solutions are given by 



U{z) = -ln[alniz)+l3] 



2 



M{z) ^ -^ (1 - 2g2) ln{[aln(z) +/?]}- ^ [aln(z) + f3]' + Mo, 

(^(z) = qln{z) + ^\n[a ln(z) + /3] + <^o- (4.12) 

2q 

From the above expressions we can see that to have the solutions valid near the hypersurface w = in the region 
M < 0, w > 0, we must assume 

a < 0. (4.13) 

To study the spacetime near the null hypersurface v = further, following the discussions given in Sec. III. B, one 
can show that now the hypersurface v — also represents a past null infinity. This null infinity is free of spacetime 
singularity. This can be seen, for example, from the Ricci scalar, which now is given by 

^ , .._ V[aln(z) + /3f~a2 



2q^zu^aHz) + pf+^''"^/^''" 
xexp|Mo-^[aln(z) + /3]H. (4.14) 

Thus, i; = is actually a past null boundary of the spacetime. On the other hand, from Eq.(4.14) we can see that 
the spacetime is singular on the hypersurface 

z = zo = e-'^/", (4.15) 

which is spacelike and serves as the up boundary of the spacetime. The corresponding Penrose diagram is given by 
Fig. 2. However, unlike the case discussed in Sec. III. B, now the scalar field is spacelike near the null hypersurface 
f = 0. Therefore, the corresponding solution cannot be interpreted as representing gravitational collapse of the scalar 
field in the region u < 0, v > 0. 



V. CONCLUSIONS 

In this paper, we have studied plane symmetric self-similar solutions to Einstein's four-dimensional theory of gravity 
and found all such solutions. We have also studied the local and global properties of those solutions and found that 
some of them can be interpreted as representing gravitational collapse of the scalar field. During the collapse, trapped 
surfaces are never developed. As a result, no black hole is formed. Although the collapse always ends with spacetime 
singularities, it is found that these singularities are spacelike and not naked. 

An interesting question is: Do these solutions represent critical collapse? To answer this question, one needs to 
study the linear perturbations of the solutions and show that there exists one solution that has one and only one 
unstable mode. Such a study is out of the scope of this paper, and we wish to return to this problem in another 
occasion. 



[1] B. J. Carr and A. A. Coley, Class. Quant. Grav. 16, R31 (1999); H. Maeda, T. Harada, H. Iguchi, and N. Okuyama, Prog. 

Theor. Phys. 108, 819 (2002); 110, 25 (2003); B. J. Carr and C. Gundlach, Phys.Rev. D67, 024035 (2003). 
[2] M. W. Choptuik, Phys. Rev. Lett. 70, 9 (1993); "Critical Behavior in Massless Scalar Field Collapse" in Approaches to 
Numerical Relativity, Proceedings of the International Workshop on Numerical Relativity, Southampton, December, 1991, 
Edited by Ray d'Inverno (Cambridge University Press, Cambridge, 1992); '^Critical Behavior in Scalar Field Collapse," in 
Deterministic Chaos in Ceneral Relativity, Edited by D. Hobill et al. (Plenum Press, New York, 1994), pp. 155-175. 
[3] C Cundlach, "Critical phenomena in gravitational collapse: Living Reviews," gr-qc/0001046 (2000), and references 

therein. 
[4] A.Z. Wang, "Critical Phenomena in Cravitational Collapse; The Studies So Far," gr-qc/0104073, Braz. J. Phys. 31, 188 

(2001). 
[5] A.Z. Wang, Phys. Rev. D68, 064006 (2003). 
[6] D. Garfinkle, Phys. Rev. D63, 044007 (2001). 

[7] F. Pretorius and M. W. Choptuik, Phys. Rev. D62, 124012 (2000). 
[8] V. Husain and M. Olivier, Class. Quantum Grav. 18, LI (2001). 
[9] D. Garfinkle and C Gundlach, Phys. Rev. D66, 044015 (2002). 
[10] E.W. Hirschmann, A.Z. Wang, and Y. Wu, "Collapse of a Scalar Field in 2+1 Gravity," gr-qc/0207121 (2002). 
[11] C R. Evans and J. S. Coleman, Phys. Rev. Lett. 72, 1782 (1994); T. Koike, T. Kara, and S. Adachi, ibid., 74, 5170 (1995); 

C. Gundlach, ibid., 75, 3214 (1995); E. W. Hirschmann and D. M. Eardley, Phys. Rev. D52, 5850 (1995). 
[12] D. Kramer, H. Stephani, E. Herlt, and M. MacCallum, Exact Solutions of Einstein's Field Equations (Cambridge University 

Press, Cambridge, England, 1980). 
[13] D. Tsoubelis and A.Z. Wang, Gen. Relativ. Grav., 22, 1091 (1990). 
[14] M.E. Cahill and A.H. Taub, Commun. Math. Phys. 21, 1 (1971). 

[15] A.H. Taub, " Plane- Symmetric Similarity Solutions for Self-Gravitating Fluid," in General Relativity, papers in honour of 
J.L. Synge, edited by L. O'Raifeartaigh (Clarendon Press, Oxford, 1972), pp. 133-150. 



